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C3 ^ We present an update of the theoretical predictions for the cross section of top-quark pair produc- 

tion at Tevatron and LHC. In particular we employ improvements due to soft gluon resummation 
at next-to-next-to-leading logarithmic accuracy. We expand the resummed results and derive an- 
alytical finite-order cross sections through next-to-next-to-leading order which are exact in all 
logarithmically enhanced terms near threshold. These results are the best present estimates for 
the top-quark pair production cross section. We investigate the scale dependence as well as the 
sensitivity on the parton luminosities. 



1 Introduction 

Top-quark pair production at the LHC is important as the collider will accumulate very high statis- 
tics for this process. In the initial low luminosity run (~ 10/(fbyear)) approximately 8-10^ top- 
quark pairs will be produced per year [1,2]. This data will allow for numerous measurements, e.g. 
of the top-quark mass, the electric charge of the top-quark or the weak couplings. Furthermore 
the data will allow precise tests of the production and the subsequent decay mechanism including 
anomalous couplings and top-quark spin correlations. See for example Ref. [3] for a recent review 
on top-quark physics at hadron colliders. 

A necessary prerequisite which all these studies share in common is, of course, a detailed un- 
derstanding of the production process. In Quantum Chromodynamics (QCD) this includes the 
radiative corrections to the cross section of heavy-quark hadro-production at the next-to-leading 
order (NLO) [4-6] together with its scale dependence as well as the dependence on the parton 
luminosities through the parton distribution functions (PDFs) of the proton. Further improvements 
of the perturbative stability through resummation of large Sudakov logarithms to next-to-leading 
logarithmic (NLL) accuracy have been considered as well [7, 8] and employed to generate ap- 
proximate results at the next-to-next-to-leading order (NNLO) in QCD [9]. Very recently also an 
estimate of bound state effects has been presented [10]. 

In particular our knowledge on the parton luminosities has constantly improved over the last 
years. Thus, it is an imminent question how these improvements in the determination of the parton 
distribution functions from global fits affect predictions for physical cross sections at LHC, which 
are sensitive to the gluon distribution function in the regime where x ~ 2.5 ■ 10^^. For LHC ob- 
servables this aspect has been quantitatively approached very recently by investigating correlations 
of rates for top-quark pair production with other cross sections [11]. For Tevatron, which to date 
has provided us with a lot of information on the top-quark, most prominently a very precise deter- 
mination [12] of its mass, m^ = 172.6 ±0.8 (stat.) ±1.1 (syst.) GeV, the cross section Opp^tix in 
Eq. ([D) had been studied some time ago [13, 14]. However, due to changes in the available PDF 
sets from global fits an update also seems to be in order here as well. 

It is the aim of this article to review theoretical predictions for the production cross sections 
of top-quark pairs at Tevatron and LHC and to establish the present theoretical uncertainty. We 
provide an update of the NLL resummed cross section as defined in Ref. [8] (and also used 
in Ref. [14]). Subsequently, we extend these results to the next-to-next-to-leading logarithmic 
(NNLL) accuracy and derive approximate NNLO cross sections thereby improving previous cal- 
culations [9, 13]. At two loops we are thus in a position to present all logarithmically enhanced 
terms near threshold and to assess their phenomenological impact by studying the quality of the 
perturbative expansion, i.e. the properties of apparent convergence and the stability under scale 
variations. This seems particularly interesting considering not only the anticipated experimen- 
tal precision at LHC [1,2] but also in view of recent activities aiming at complete NNLO QCD 
predictions for heavy-quark hadro-production [15-19]. 

The paper is organized as follows. In Sec. [21 we set the stage and study the threshold sensitivity 
of the inclusive hadronic cross section for top-quark pair production. Subsequently, we provide 
updates of Refs. [8, 14] employing recent sets of PDFs. In Sec. [3] we extend the resummed cross 
section to NNLL accuracy and calculate the complete logarithmic dependence of the cross section 
near threshold (including the Coulomb corrections). Together with the exact NNLO scale depen- 
dence of Ref. [9] these results are the best present estimates for the hadro-production cross section 



of top-quark pairs. We conclude in Sec.|4]and give some relevant formulae to the Appendix. In ad- 
dition cross sections predictions using different approximations for individual PDFs are also listed 
in the Appendix. 

2 Theory status 

Throughout this article, we restrict ourselves to the inclusive hadronic cross section Opp^tix (see 
e.g. Ref. [20] for recent work on top-quark pair invariant mass distributions). Denoting the 
hadronic center-of-mass energy squared by 5had and the top-quark mass by m^ the total hadronic 
cross section for top-quark pair production is obtained through 
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The parton luminosities L/v(f,5had7A'f) are defined through 
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where fi/p{x,^j) is the PDF describing the density of partons of flavor / in the proton p carrying a 
fraction x of the initial proton momentum, at factorization scale /jr. Note that we have included .had 
into the definition of Ljj to allow an easy comparison of the luminosity function between different 
colliders, in particular LHC and Tevatron. The sum in Eq. ([T]) runs over all massless parton flavors 
and the top-quark mass used is the so-called pole mass. 

Within the context of perturbative QCD the standard way to estimate the theoretical uncertainty 
for the inclusive hadronic cross section Opp^ttx in Eq. ([T]) is based on the residual dependence on 
the factorization/renormalization scale /j^//j^. Starting from the available predictions to a certain 
order in perturbation theory it is common practice to identify the factorization scale with the renor- 
malization scale (i.e. /jr = fi^ = /u) and to estimate the effect of uncalculated higher orders by 
varying ^ in the interval [wf /2, 2mt] . For a given global PDF fit in contrast, the uncertainties which 
stem from uncertainties of the experimental data used in the fits are treated systematically by a 
family of npop pairs of PDFs, where npoF is the number of parameters used in the fit. Then, the 
systematic uncertainty for the observable O under consideration is estimated by (e.g. [11,21]), 

^0^1 J E (0,+ -0,_)2. (3) 
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Here the observable Ok± are obtained by using the parton distribution functions f., ^ obtained 
by a "statistical" ± la- variation of the fcth fit parameter after diagonalization of the correlation 
matrix. (Strictly speaking the fit parameters are varied by an amount which the authors of the 
corresponding PDF set take to be equivalent to a ± la- variation.) To end up with an estimate 
of the overall uncertainty the uncertainty coming from the PDFs has to be combined with the 
uncertainty due to uncalculated higher orders. Given that the two uncertainties are very different 



from each other: in one case we are faced with the traces of an experimental uncertainty — in the 
other case we have the systematic uncertainty due to missing higher order corrections which clearly 
do not follow any statistical law. Adding in quadrature the two uncertainties therefore seems to be 
inappropriate and we use a linear combination of the uncertainties as a conservative estimate for 
the total uncertainty of the top-quark pair cross-section: 

a{2mt) -A(5pDF{2mt) < o < o(mf/2) +AopDF{mt/2) , (4) 

where AopoF is computed according to Eq. ([3]). The NLO QCD corrections for the partonic cross 
sections Otj^tt in Eq. Qi (known since long time [4-6]) provide the first instance in this procedure 
where a meaningful error can be defined through Eq. ©. 

Thus, let us start our discussion of the various theoretical uncertainties by reviewing at NLO 
in QCD some basic aspects concerning the parton luminosities Ltj as defined in Eq. (O at LHC 
(Fig.d]) and Tevatron (Fig.O. At LHC the highest flux is provided by the quark-gluon initial state 
(first plot of Fig. [T]). However, as can be seen from the cross section plot of Fig. [H the parton level 
cross section Oqg^tt is much smaller than for the qq- or ^^-initiated processes, since the ^g-channel 
is of order a^-^ and thus formally a NLO correction. As a consequence, in the total hadronic cross 
section the (^g-channel gives only a contribution at the percent level (see second last plot in Fig. [T]). 
In principle the same argument applies for (^g-channel, however its contribution to the hadronic 
cross section is even further suppressed because of the smaller parton luminosity. 

The second largest parton flux is delivered by the gg-channel. In combination with the large 
partonic cross section Ogg^tt this is the most important channel at LHC, resulting in about 90% 
of all top-quark pairs produced via gluon fusion. Close to threshold the uncertainty of the gluon 
flux (estimated using Eq.[3l see also discussion there) is about 3% and at 1 TeV it grows to almost 
10%. However, given that the gg-channel is largely saturated at parton energies v^ ~ 1 TeV, the 
large PDF uncertainty above 1 TeV does not have significant impact on the overall uncertainty of 
^pp->ttx in Eq. ([T]), see last plot in Fig.[TJ 

The parton luminosity Lqq ranks third at LHC. Close to threshold the qq-^ux is suppressed by 
roughly a factor 10 compared to Lgg. The corresponding parton cross section dqq^,j vanishes in 
the high energy limit in contrast to the gg- and ^g-case where the partonic cross sections approach 
a constant at high energies. Thus, the ^^-contribution to the hadronic cross section saturates well 
below 1 TeV and adds the known 10% at LHC. Due to the small numerical contribution of the qq- 
and ^g-channels the PDF uncertainty of the hadronic cross section Opp^ttx is entirely dominated 
by the uncertainty of Lgg. 

At Tevatron, the situation is reversed (see Fig.|2l). The luminosities Lij are ordered in magnitude 
according to Lqq > Lqg > Lgg. This makes the ^<^-channel by far the dominant one contributing 85% 
to the hadronic cross section Opp^tfx, while gluon fusion almost makes up for the rest. Although 
the PDF uncertainty of the qq-Qux are only 3-4% at low energies the overall PDF uncertainty of 
the top-quark cross section Opp^ttx is large, because of the sensitivity to the gluon PDF content at 
large-x, which is still poorly constrained at present (see Fig. |2l). 

Finally, it is interesting to determine the value of Smax for which the cross section 
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Figure 1 : The parton luminosity L,y with the individual PDF uncertainties (upper plots) and the parton 
cross sections dtj^tt at NLO in QCD (third plot from below) as a function of the parton energy ^/s. The 
lower plots scan the total cross section Cj(i'had,'Mf^;i'max) (the total PDF uncertainty) as a function of y/Smax 
for LHC. We use .Js^ = 14 TeV, m, = 171 GeV, ^ = m^ and the CTEQ6.5 PDF set. The dashed Une 



indicates the value of ^/s^^ for which the cross section is saturated to 95%. 
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Figure 2: Same as Fig. [T]for Tevatron {-^/s^a = 1 -96 TeV). 



saturates the total cross section a(5had, "^f ^) to 95%. At Tevatron this happens at y/Smax ~ 600 GeV 
as can be seen from Fig. [2l Thus, the total cross section is largely dominated by parton kinematics 
in the range -\/f ~ 2mj close to the threshold of the top-quark pair. This makes top-quark pair 
production at Tevatron an ideal place to apply threshold resummation. At LHC energies in contrast, 
the available phase space is larger and saturation to 95% is only reached at parton energies ^5max ~ 
1 TeV (see Fig. \Q. This makes the cross section less sensitive to Sudakov logarithms, although 
numerically a significant part still originates from the threshold region for parton kinematics due 
to the steeply decreasing parton fluxes. 

Further refinements of perturbative predictions for Opp^ttx in Eq. ([T]) do rely on subsequent 
higher orders to be calculated. In particular, the knowledge about large logarithmic corrections 
from regions of phase space near partonic threshold allows for improvements of the theoretical 
accuracy beyond NLO in QCD. These Sudakov type corrections can be organized to all orders 
by means of a threshold resummation (e.g. to NLL accuracy [7, 8]), which has been the basis for 
phenomenological predictions employing a resummed cross section as defined in Ref . [8] (and also 
usedinRef. [14]). 

However, before updating the NLL resummed results of Ref. [8] let us briefly give some rele- 
vant resummation formulae. It is well known that soft gluon resummation for ?f-production relies 
on a decomposition of the parton-level total cross-section in the color basis, conveniently defined 
by color-singlet and color-octet final states. Then we can decompose 

dij^tt{s,mt^,fjj,fj^)= Y, ^ij,i{s,mt^,Mf,f^r)- (6) 

7=1,8 

Moreover, we use the standard definition of Mellin moments 

1 

dfjj{mt^,fi},p^) = J dpp^-^6ijj{p,mt^,^ij,^i^), (7) 







with 



P = ^. (8) 



Then, the resummed Mellin-space cross sections (defined in the MS-scheme) for the individual 
color structures of the scattering process are given by a single exponential (see e.g. Refs. [22,23]), 



%'!^' :'':''] = 4,/K''A/?,A/.^) -exp (G^Hm.',plp^)) +0{N-^ ln"N) , (9) 

where 6j: } denotes the Bom term and the exponents Gf: j are commonly expressed as 

Gfjj = \nN-glj{l)+gl^,{l)+asglj{l) + ..., (10) 

where X = ^oaslnN and as — as/{4n). To NLL accuracy the (universal) functions g]- as well 
as the functions g? j are relevant in Eq. (flOl) . of course, together with the appropriate matching 
functions g^j j in Eq. Q. Explicit expressions can be found below. 



For phenomenological applications [8] the soft-gluon resummation in A'^-space at the parton 
level one introduces an improved (resummed) cross section Ores, which is obtained by an inverse 
Mellin transformation as follows, 
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where ^^l^jj is the standard fixed order cross section at NLO in QCD and d^ /Inlo ^^ ^^^ perturba- 
tive truncation of Eq. ^ at the same order in a.;. Thus, the right-hand side of Eq. (fTTl) reproduces 
the fixed order results and resums soft-gluon effects beyond NLO to NLL accuracy. 
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Figure 3 : The partonic cross sections for the processes qq — > tt (left) and gg -^ tt (right) in pb, and for 
fj = Mt = 171 GeV. The dotted lines ai^e the exact NLO result [4, 5] and the solid lines correspond to the 
NLL resummed cross sections [8] (see text for details). 



In Fig. [3] we plot the resummed cross section Ores//^/F as defined in Eq. (fTTl) zooming in on 
the threshold region. We display the qq-channe\ (left) and the gg-channel (right) as a function of 
the distance from the partonic threshold at v^= Init. As mentioned already the qg-channel is sup- 
pressed by an additional power of a.^ and therefore does not contribute large Sudakov logarithms 
to the accuracy considered here. 

The results for Ores as shown in Fig. [3] have been obtained by performing the inverse Mellin 
transform in Eq. (fTTl) numerically. Following the procedure described in Ref. [8] we have com- 
pared our results obtained from the numerical inversion with the ones shown in Ref. [8] and found 
complete agreement. To be precise, the treatment of the constant terms in Ref. [8] (i.e. the terms 
denoted g^- j in Eq. (|9|)) differs slightly from the minimal approach, e.g. Eq. (fTTl) . Some constants 
which are formally subleading have been included in Ref. [8], and moreover, several schemes for 
power suppressed terms in A'^ have been implemented. The resummed result shown in Fig. |3] is 
defined through Eq. (63) of Ref. [8] with the parameter A set to 2. Another issue concerns the 
precise numerical matching of the exact NLO cross section and the resummed result in Eq. (fTTl) . 
We apply the resummed result only for \f§—2mt < 10 GeV. The exact point is determined from 



the crossing of the two curves for Ores and Ojv^lq. Note that the precise numerical value is not 
important. Effectively Ores is thereby restricted to a region of parton energies of V^ ~ 2mt or in 
other words, to a kinetic energy of one top-quark of a few GeV. 
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Table 1 : The NLL resummed cross section of Ref. [8] in pb (see text for details) for various values of 
the top-quark mass m, at Tevatron (-^^had = 1-96 TeV) using the CTEQ6.5 PDF set [24]. 6 is the relative 
uncertainty with respect to the central value: 5 = 100 x (max — min) /(max + min). 



It should be stressed, that the total hadronic cross section is not very sensitive to these fine 
details due to the convolution with the parton luminosities L/y in Eq. ([U). However, to per mille 
accuracy (see Tabs.[l]-l4l) these details become noticeable. Having them clarified, we are now in a 
position to update previous results [8, 14] for the rf-cross section using modem PDFs, such as the 
CTEQ6.5 PDF set [24]. In comparison to older sets, we find e.g. a shift of 3% in the total cross 
section between the PDF sets CTEQ6.5 [24] and CTEQ6.1 [21], see also Ref. [11]. Since Eq. (dU) 
effectively contains the dominant part of higher orders (NNLO and beyond) it seems however 
equally appropriate to use also the MRST-2006 NNLO PDF set [25]. From the results in Tabs.[I]-l4] 
and Fig. m we conclude that the present overall uncertainty on the NLL resummed ff-cross section 
is 12% at Tevatron with a small shift of the central value between different PDF sets. At LHC 
the NLL resummed cross section of Ref. [8] reduces effectively to the NLO QCD prediction of 
Refs. [4, 5] with an overall uncertainty of 14%. This is due to the small gluon contribution close 
to threshold (see Fig. [3]) and the cut when matching NLO and the Ores in the numerical determi- 
nations. Thus the reduction of the theoretical uncertainty through Eq. (fTTl) is marginal in this case. 
We will see in the next Section [3l how this situation can be improved with the help of approximate 
NNLO QCD corrections. 
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Table 2: Same as in Tab. fusing the MRST-2006 NNLO PDF set [25]. 
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Table 3: The NLL resummed cross section of Ref. [8] in pb (see text for details) for vaiious values of the 
top-quark mass nit at LHC (Js^ = 14 TeV) using the CTEQ6.5 PDF set [24]. 
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Table 4: Same as in Tab. fusing the MRST-2006 NNLO PDF set [25]. 

3 Prospects at NNLO in QCD 

Let us now extend the theory predictions for heavy-quark hadro-production. We will focus on the 
threshold region and improve soft gluon resummation to NNLL accuracy. Subsequently, we em- 
ploy the resummed cross section to generate higher order perturbative corrections - more specifi- 
cally an approximate NNLO cross section CJjyfjyjLQ (apnrox) which is exact to logarithmic accuracy 
(including the Coulomb corrections). To that end, we briefly recall the steps leading to the final 
form for G^- j in Eq. (flOl) . In order to achieve NNLL accuracy the function g? j is of particular 
interest here. 

The exponential G^ j in Eq. ^ is build up from universal radiative factors for the individual 
color structures which take the form 






^DY/Higgs + ^A 






(12) 



where the exponentiation of singlet contribution (i.e. a colorless massive final state) follows from 
the Drell-Yan process and hadronic Higgs production in gluon fusion. The corresponding functions 
Gqy and G^j are very well-known [26-28]. The exponentiation of the color-octet contribution 
receives an additional contribution G^^ due to soft gluon emission from the heavy-quark pair in 

the final state. The final-state system carries a total color charge given by the QQ charge, thus its 
contribution to soft radiation vanishes in the color-singlet channels regardless of the initial state 
partons. IVIoreover gluon emission from massive quarks does not lead to coUinear logarithms and 
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Figure 4: The tt total cross section resummed to NLL accuracy [8] as a function of nit for the Tevatron at 
Y^^had = 1-96 TeV (left) and LHC at y^^had = 14 TeV (right). The solid line is the central value for /u = nit, 
the dashed lower and upper lines coiTcspond to ^ = 2m, and ^ = m,/2, respectively. The band denotes the 
total uncertainty that is the uncertainty due to scale variations and the PDF uncertainty of the CTEQ6.5 
set [24] combined together according to Eq. @. 



therefore G^^ starts at NLL accuracy only. Explicit formulae are 



^DY/Higgs 



r 7^-1 - 1 r4m,2(l-z)2 J„2 



',} 



/7 — 1 

dz^-j-^DQQiasiW[l-z]^)) 



(14) 



with anomalous dimensions A, and D,, i = q,g corresponding to DY and Higgs, respectively. The 
effects of coUinear soft-gluon radiation off initial-state partons i = q,g are collected by the first term 
in Eq. (fT3l) while the process-dependent contributions from large-angle soft gluons are resummed 
by the second term. Soft radiation from the heavy-quark pair in the final state is summarized by 
in Eq. (fT4l) with the corresponding anomalous dimension Dqq. 
The extension to NNLL requires the Aq and Ag to three loops [29, 30] and the function Dq and 
Dg to two loops [26-28] (the latter are actually known to three loops as well [31,32]). Explicit 
expressions using the expansions 






fiCLs 



y f{i) ^ 

V^ 4:r 



1/ 



«4. 



(15) 



are collected in Eqs. (iA[2l) - (IA[4l) in the Appendix. The remaining anomalous dimensions Dqq for 
soft radiation off a heavy-quark pair in the final state is needed to two loops. For the latter we use 
the exact calculation of the two-loop QCD corrections to the massive heavy quark form factor [33] 
along with Ref. [34], where the exponentiation of the form factor for massive colored particles in 
the limit m^ -^ has been clarified. From the all-order singularity structure of the massive form 
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factor [34], we can read off the single poles corresponding to soft gluon emission at one and two 
loops. After a trivial substitution of color factors (Ca for Cf) we find 



D 



QQ 



D 



,(2) 



QQ 



(16) 



As pointed out in Ref. [8] the one-loop value of Dqq agrees with Ref. [7] where the soft anomalous 
dimension matrix (in color space) Tjj for heavy-quark production has been calculated at order a^. 
In the limit (3 — > the matrix r)^ diagonalizes in the singlet-octet basis and reproduces D j^ from 
its eigenvalue in the octet channel. Moreover, the structure of Dqq as determined from Ref. [34] 

(2) 

agrees also with the two-loop soft anomalous dimension matrix r)y for vanishing parton masses 
which obeys the following factorization property [35, 36], 



IJ 



m=0 



Usl ij 



m=0 



^(2) 



(17) 



with the well-known ratio Ag /Ag [37]. As a further check on Eq. ([T6l) it would, of course, be 
very interesting to repeat the calculation of Ref. [35, 36] for heavy-quark hadro-production at two 
loops, i.e. with non-vanishing parton masses m ^ 0. 
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Figure 5: The resummed cross section in Mellin space for the qq (left) and the gg-channel (right) normal- 
ized to the Born result for /j = wi, = 171 GeV. The dotted lines are the LL approximation, the dashed lines 
denote the NLL result and the solid lines coiTespond to the NNLL result derived in this paper. 



Finally, explicit integration of Eqs. (fT3])-(IT4l) leads to the functions gj-, gj- j, g^- j of Eq. ([10]) 
and to the matching g^- j in Eq. ^ which we collected in Eqs. (IA[5I) - (IA[13I) in the Appendix. 
All formulae can be obtained by simple substitutions e.g. from deep-inelastic scattering (DIS) in 
Ref. [28]. In Fig. [5] we display the resummed cross section of Eq. ^ (normalized to the respective 
Born result) for the qq and the gg-channel in increasing logarithmic accuracy. Fig. |5] clearly shows 
the good convergence property of the NNLL contribution similar to other observables investigated 
previously [26-28]. As a matter of fact, Eq. (flOl) may even be extended to next-to-next-to-next- 



to-leading accuracy (N-^LL), as the relevant functions g^ ^ can be easily derived from Ref. [28] 
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and the respective anomalous dimensions are known. Following the arguments based on the ex- 
ponentiation of the form factor for massive colored partons [34] leading to Eq. (fT6l) we identify 

D^ = —Ag and for the four- loop terms Aq and Ag a Pade estimate exists [28]. However, 
presently we are lacking knowledge on the matching functions g^: at this order. From experience 
we expect N^LL effects to be numerically very small, though, and we leave this issue to future 
investigations. 

Let us instead use the resummed cross section Ores (now known at NNLL accuracy) to con- 
struct an approximate NNLO cross section ctjvjjvjlq (apnrox) ^^ expanding Eq. ^ to second order. 
In this way, we determine the corresponding Sudakov logarithms appearing in the NNLO correc- 
tions, i.e. the powers of In^A'^ in Mellin space or In^ (3 in momentum space with k = 1 , ..., 4 and the 
velocity of the heavy quark 



(3=y^l-4m,2/^. 

As a technical remark, we stress here that in a fixed order expansion the inverse Mellin transfor- 
mation, i.e. the mapping of powers of InA'' to powers of ln((3) can be uniquely performed. The 
Mellin-space accuracy up to power suppressed terms in A'^ corresponds to neglecting higher order 
polynomials in (1 — p) with p = 4mt^/s. We give some formulae in the Appendix. Moreover, we 
can even include the complete Coulomb corrections at two loops thanks to Ref. [33,38]. In the 
singlet-octet decomposition dijj of the individual color structures in the cross section, the result 
of Ref. [33, 38] can be directly applied to the singlet case, while a simple modification of the color 
factors (that is (C/r — Ca/2) instead of C^) accounts for the octet case (see e.g. [39]). 

Thus, we are in a position to present the threshold expansion of the inclusive partonic cross 
sections Gij^tt{s,m^,i-i^) entering Eq. ([I]). In the perturbative expansions in powers of the strong 
coupling constant a^ as defined in Eq. ([T5l) and setting /u = nit and nf = 5, we have for the qq 
channel in the MS -scheme 



-(1) 
o - - 

qq^tt 



-(2) 
O - - 
qq^ttt 



aj,^) /42 667in2p-20.6101n(3+13.910-3.2899H, (18) 

qq Hi I [^ I 

^SL/f| 910.221n4 p - 1315.51n3p+ (565.80- 140.37^ j In^p 



.^j In p + 3.6077^+10.474^ ' ^^^^ 



+ (862.42 + 32.106- ) In p + 3. 6077-2 + 10.474- +qj^ [> , (19) 



o^P = a[,"j_^„-|961n2p-9.51651np + 35.322 + 5.1698^J>, (20) 



o 



gg^tt gg->tt 

(2) 



^tt 



a^o)^^-|46081n4|3- 1894.91n3(3+f-3.4811 +496.30^ jln^p 

p; " p^ ■ p 



+ ( 3144.4 + 321.17- ) In |3 + 68.547^- 196.93- +Qg \ , (21) 



The terms proportional to inverse powers of P correspond to the Coulomb corrections and the 

(2) (2) 

presently unknown two-loop constants Q= and Cgg are set to zero. For reference, we have also 
repeated the well-known NLO results Eq. ([TSl) and Eq. (I20l) . The lengthy analytical results (con- 
taining the explicit dependence on the color factors Ca^Cp and on Uf) are given in the Appendix, 
Eqs. (IA[17I) - (IA[20I) . We define a NNLO (approx) cross section to be used in this paper as the sum 
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of exact NLO result and the two-loop contribution of Eqs. ( 1191 ) and (1211) for all scale independent 
terms. All scale dependent terms at NNLO accuracy are long known exactly [9] as they can be 
easily constructed from the lower orders convoluted with the appropriate splitting functions. We 
use the exact result of Ref. [9] for the /j-dependence at two loops. As emphasized several times, the 
qg- and ^^-contributions are small at Tevatron and LHC and we simply keep them at NLO here. 
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Figure 6: The partonic cross sections for the processes qq — > tt (left) and gg — > tt (right) in pb, and for 
/Li = Mt = 111 GeV. The dotted lines ai^e the exact NLO result [4, 5] and the solid lines coiTcspond to the 
NNLO (approx) result of this paper. The dashed lines close by denote the previous approximations of [9]. 



In Fig. [6] we display the NNLO (approx) results for the partonic cross sections as calculated 
in Eqs. (fT9l) and (f2T]) . We plot Gqq-,tt (left) and Ogg^tt (right) again as a function of the distance 
from threshold at a/? = 2mt (solid lines in Fig. (6]). This is useful to asses the question what 
improvements can be expected from a full NNLO calculation. Given that the ^-factor of the NLO 
result is of moderate size at large partonic energies one can expect that the full NNLO corrections 
should give only small corrections in this region of phase space. On the other hand the ^-factor 
of the NLO correction becomes large in the threshold region (i.e. \f§ — 2mt < (50) GeV), where 
finally perturbation theory breaks down. In this region the corrections are dominated by the large 
logarithms in P together with the Coulomb corrections which are correctly described by Eqs. ( fT9l ) 
and (|2TI) . In other words, where we expect large corrections from the full NNLO QCD calculation 
our CTjvjjvjLQ (approx) should provide a good estimate. This is further supported by Ref. [15] where 
the next-to-leading order corrections to top-quark pair production together with an additional jet 
have been calculated. This contribution represents part of the NNLO corrections for inclusive top- 
quark pair production. In Ref. [15] it was found that for ^ = ^r = ^^ = mt the NLO corrections to 
tt+ 1-jet production are almost zero. This is a further indication that the hard corrections to the 
inclusive top-quark pair production at NNLO are indeed small. IMoreover, as mentioned above, we 
have further improved ctjv^js^lq (approx) ^y incorporating the complete scale dependence at NNLO 
which is already known exactly [9]. To that end let us quantify in detail once more the range of 
validity for the soft gluon approximations. We show in Tab. [5] the numerical size of the individual 
logarithms and powers in (3 in Eqs. (fT8l)-(l2TI) when evaluated at a given parton energy s. As a 
consequence of the moderate size of the expansion coefficients in Eqs. (fT8l)-(l2T]). we clearly see 
the good convergence properties of the logarithmic expansion up energies \/§ — 2mt < (50) GeV. 
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Beyond that value the soft logarithms (being proportional to the Born cross section) smoothly 
vanish. 
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50.0 


0.2628 


-0.3671 


0.5127 


-0.7160 


4.1870 


2.0462 


0.4887 
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1.0943 
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Table 5: Numerical values of the individual powers of InP and P for various distances from threshold 
\/l—2nif for rrit = 171 GeV as entering in Eqs. (fT8])-(l2TI). 



In addition we show in Fig. [6] also previous approximations to the NNLO correction from 
Ref. [9] (dashed lines) employing two distinct differential kinematics to define the partonic thresh- 
old. They agree with our NNLO (approx) corrections for \f§—2mt < 30 GeV if Eqs. (fT9l) and (|2TI) 
are truncated to the first three powers in ln|3. However, at higher partonic center-of-mass energies 
the results of Ref. [9] receive large numerical contributions from sub-leading terms and become 
unreliable (see in particular Fig. [6] on the right). 

We would also like to point out that there is a discrepancy between the NLL resummed cross 
section of Ref. [8] and Eq. (fTTI) and the fixed order NNLO approximation discussed here. In par- 
ticular for the gluon fusion channel in Fig. [3] (right) and Fig. [6] (right) the numerical differences 
between Ores and o^NLO (approx) ^^^ rather large. The all-order NLL resummed cross section 
Ores is significantly smaller than its expansion to second-order Eq. (f2T|) or the corresponding result 
of Ref. [9], the latter two both being consistent with each other. We can attribute this difference to 
the following fact: For rf-hadro-production the Born cross section exhibits simple (although non- 
trivial) A/^-dependence and the resummed cross section Ores in Mellin space (as implemented in 
Ref. [8] and Eq. ([TT]) ) contains products of A'^-dependent functions. This is unlike other cases con- 
sidered in the literature, where the Born terms have always been proportional to a delta-function, 
i.e. 5(1 —x) for DIS, Drell-Yan or Higgs production. In momentum space products of A'^-dependent 
functions correspond to convolutions, which induce formally sub-leading but numerically large 
corrections in the resummed result. Eventually, this leads to the observed suppression in Fig. [3l 
Most likely this large discrepancy will also persist when comparing Ores to a full NNLO QCD 
calculation, or upon matching the latter to a resummed cross section at NNLL accuracy along the 
lines of Eq. (fTT)) . This fact has to be kept in mind when using Ores for predictions at LHC, where 
the gg-channel dominates. As mentioned above, for any finite order expansion in a.?, there is no 
ambiguity in performing the inverse Mellin transformation analytically up to power suppressed 
terms in A^ or, equivalently in ( 1 — p) . 

We are now in a position to present the new results for the top-quark cross section Opp^ttx at 
NNLO (approx) as defined below Eq. (|2TI) including the exact scale dependence. We also quote 
the corresponding uncertainty according to Eq. dU. In our study we use the same PDFs as in 
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Figure 7: Comparison of the scale dependence of c^nnLO (approx) ^^^^ ^^^ ^^^ MRST-2006 NNLO [25] 
and CTj«^LO ^^^^ ^^^ ^^^ CTEQ6.5 [24]. For comparison we show also the corresponding results for the 
Alekhin set of PDFs [40]. The cross sections are normalized to the value at/^ = m^. 



Tabs. [T]-l4] above. In Figj7]the scale dependence for Oj^lq and Oj^js^lq (approx) ^^ shown. For the 
^NNLO (approx) we use the PDF set ]y[RST-2006 NNLO while for Onlo *^ ^^^ ^et CTEQ6.5 
is used. To become less sensitive to different normalizations we normalise the curves to the central 
value (5{iJ = Mt). In addition we show also the results obtained by using the Alekhin PDF set [40] 
available in NLO and NNLO accuracy. After normalisation one can see that the two curves for the 
NLO predictions agree rather well as one might expect. For the two NNLO curves the agreement 
is less good in particular for extreme values of the scale. The origin of the minor discrepancy 
might be attributed to slightly different input densities at low scale. Compared to the NLO results 
the scale dependence of cJj^nLO (approx) ^^ improved. In particular we find a plateau at ^u = m^ . 
Note that due to the different shape of the NNLO curve the uncertainty estimate of Eq. Q has to 
be adapted here. Restricting the scale to the interval [mt/l.lnit] the residual scale dependence of 
^NNLO (approx) ^^ reduced to a few percent. 

In Fig. [8] we plot the total cross section at Tevatron and display also CDF data [41] with 
mj = 171 GeV. The corresponding cross section values are given in Tabs. [6] and |7l By using 
^NNLO (approx) ^^^ residual scale dependence is reduced to 3%. Compared to Ojs^lq and Ores 
presented in the previous section this corresponds to a reduction by a factor of 2. The data points 
nicely agree with the theoretical prediction. The overall uncertainty of the the theoretical pre- 
diction is about 8% for CTEQ6.5 and 6% for lMRST-2006 NNLO — an accuracy unlikely to be 
reached at the Tevatron experiments. Note that the smaller PDF uncertainty obtained when using 
the lMRST-2006 NNLO PDF set is due to a different convention used by the IMRST collaboration 
to define the PDF uncertainty. In Fig. |9] (right) we show ctnNLO (approx) ^°^ ^^^ LHC. Again we 
observe a drastic reduction of the scale uncertainty compared to cjj«^lq and Ores- For comparison 
we show also ctnLO ^^ ^^S- l2l(left)- The corresponding cross section values are listed in Tabs. [8] 
and[9l Apart from reducing the scale dependence the o^NLO (approx) le^ds only to a small shift 
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of a few percent in the central value for the cross section prediction. The Oj^f^LO (approx) band 
is well contained in the Ojs^lq band. So perturbation theory seems to be well behaved and un- 
der control. The overall uncertainty is about 6% for the CTEQ6.5 PDF set and about 4% for the 
MRST-2006 NNLO set. 

The numbers quoted in Tabs. |6]-l9] represent presently the best estimates for the top-quark pro- 
duction cross section at Tevatron and LHC (see the Appendix for additional information on the 
individual PDFs and their eigenvalues). It should be kept in mind, though, that there is an intrinsic 
uncertainty in the central value at /j = nit of our NNLO (approx) result due to neglected power 
corrections in |3 ~ (1 — p) away from threshold. However, due to the steeply falling parton flux 
(see Figs. [HO, the numerical impact of these contributions is much suppressed. 
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Figure 8: The NNLO (approx) QCD prediction for the ?? total cross section at Tevatron and CDF data [41] 
with m, = 171 GeV - as functions of m, for ^/si^ = 1.96 TeV (left) and of y^ishad (right). The solid line is the 
central value for ^A = mt, the dashed lower and upper lines coiTcspond to /j = 2m, and n = rnt/2, respectively. 
The band denotes the total uncertainty that is the uncertainty due to scale variations and the PDF uncertainty 
of the MRST-2006 NNLO set [25] combined together according to Eq. dH). 
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8 


174 


6.97 


7.3 


3 


6.83 


7.63 


6 


6.59 


7.69 


8 


175 


6.76 


7.08 


3 


6.62 


7.4 


6 


6.39 


7.45 


8 


176 


6.55 


6.87 


3 


6.43 


7.17 


6 


6.2 


7.23 


8 


177 


6.36 


6.66 


3 


6.23 


6.96 


6 


6.01 


7.01 


8 


178 


6.16 


6.46 


3 


6.05 


6.75 


6 


5.83 


6.8 


8 


179 


5.98 


6.27 


3 


5.87 


6.54 


6 


5.66 


6.6 


8 


180 


5.8 


6.08 


3 


5.69 


6.35 


6 


5.49 


6.4 


8 



Table 6: The cross section c^nnlq (annrox) ^^ derived in this paper in pb for various values of the top- 
quark mass nit at Tevatron (^/5'had = 1-96 TeV) using the CTEQ6.5 PDF set [24]. 





only ! 


5cale uncertainty 


only pdf uncertainty 


total 


uncertainty 


m 


min 


max 


5[%] 


min 


max 


5[%] 


min 


max 


5[%] 


165 


9.59 


10 


3 


9.73 


10.2 


3 


9.34 


10.3 


6 


166 


9.28 


9.76 


3 


9.42 


9.94 


3 


9.04 


10 


6 


167 


8.99 


9.45 


3 


9.12 


9.62 


3 


8.75 


9.7 


6 


168 


8.7 


9.16 


3 


8.83 


9.31 


3 


8.47 


9.39 


6 


169 


8.42 


8.87 


3 


8.55 


9.02 


3 


8.2 


9.1 


6 


170 


8.16 


8.59 


3 


8.28 


8.73 


3 


7.94 


8.81 


6 


171 


7.9 


8.32 


3 


8.02 


8.46 


3 


7.69 


8.53 


6 


172 


7.65 


8.06 


3 


7.77 


8.19 


3 


7.45 


8.27 


6 


173 


7.41 


7.81 


3 


7.53 


7.93 


3 


7.22 


8.01 


6 


174 


7.18 


7.57 


3 


7.29 


7.69 


3 


6.99 


7.76 


6 


175 


6.95 


7.34 


3 


7.07 


7.45 


3 


6.77 


7.52 


6 


176 


6.74 


7.11 


3 


6.85 


7.22 


3 


6.57 


7.29 


6 


177 


6.53 


6.89 


3 


6.64 


6.99 


3 


6.36 


7.07 


6 


178 


6.33 


6.68 


3 


6.43 


6.78 


3 


6.17 


6.85 


6 


179 


6.13 


6.48 


3 


6.24 


6.57 


3 


5.98 


6.64 


6 


180 


5.95 


6.28 


3 


6.05 


6.37 


3 


5.8 


6.44 


6 



Table 7: Same as in Tab. fusing the MRST-2006 NNLO PDF set [25]. 
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only scale uncertainty 


only pdf uncertainty 


total 


uncertainty 


m 


min 


max 


5[%] 


min 


max 


5[%] 


min 


max 


5[%; 


165 


1035 


1082 


3 


1048 


1117 


4 


1003 


1117 


6 


166 


1007 


1052 


3 


1019 


1086 


4 


975 


1086 


6 


167 


979 


1024 


3 


991 


1056 


4 


948 


1056 


6 


168 


953 


996 


3 


964 


1028 


4 


922 


1028 


6 


169 


927 


969 


3 


937 


1000 


4 


897 


1000 


6 


170 


902 


943 


3 


912 


973 


4 


873 


973 


6 


171 


878 


917 


3 


887 


947 


4 


849 


947 


6 


172 


855 


893 


3 


863 


922 


4 


827 


922 


6 


173 


832 


869 


3 


840 


898 


4 


805 


898 


6 


174 


810 


846 


3 


818 


874 


4 


783 


874 


6 


175 


789 


824 


3 


796 


851 


4 


762 


851 


6 


176 


768 


802 


3 


775 


829 


4 


742 


829 


6 


177 


748 


781 


3 


755 


808 


4 


723 


808 


6 


178 


729 


761 


3 


735 


787 


4 


704 


787 


6 


179 


710 


741 


3 


716 


767 


4 


686 


767 


6 


180 


692 


722 


3 


698 


747 


4 


668 


747 


6 



Table 8: The cross section c^nnlq (annrox) ^^ derived in this paper in pb for various values of the top- 
quark mass nit at LHC {^/s^ = 14 TeV) using the CTEQ6.5 PDF set [24]. 





only scale uncertainty 


only pdf uncertainty 


total 


uncertainty 


m 


min 


max 


5[%] 


min 


max 


5[%] 


min 


max 


5[%; 


165 


1094 


1141 


3 


1128 


1154 


2 


1082 


1154 


4 


166 


1064 


1110 


3 


1097 


1122 


2 


1052 


1122 


4 


167 


1035 


1080 


3 


1067 


1092 


2 


1024 


1092 


4 


168 


1008 


1050 


3 


1038 


1063 


2 


996 


1063 


4 


169 


981 


1022 


3 


1010 


1034 


2 


969 


1034 


4 


170 


955 


995 


3 


983 


1007 


2 


943 


1007 


4 


171 


929 


969 


3 


957 


980 


2 


918 


980 


4 


172 


905 


943 


3 


932 


954 


2 


894 


954 


4 


173 


881 


918 


3 


907 


929 


2 


871 


929 


4 


174 


858 


894 


3 


883 


905 


2 


848 


905 


4 


175 


836 


871 


3 


860 


882 


2 


826 


882 


4 


176 


814 


848 


3 


838 


859 


2 


804 


859 


4 


177 


793 


826 


3 


816 


837 


2 


783 


837 


4 


178 


773 


805 


3 


795 


815 


2 


763 


815 


4 


179 


753 


785 


3 


775 


795 


2 


744 


795 


4 


180 


734 


765 


3 


755 


774 


2 


725 


774 


4 



Table 9: Same as in Tab. [fusing the MRST-2006 NNLO PDF set [25]. 
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Figure 9: The NNLO (approx) QCD prediction for the tt total cross section at LHC as functions of m, 
for ^/Shad = 14 TcV (right). The solid line is the central value for /u = mt, the dashed lower and upper 
lines correspond to /i = Inif and /j = wi,/2, respectively. The band denotes the total uncertainty that is 
the uncertainty due to scale variations and the PDF uncertainty of the MRST-2006 NNLO set [25]. For 
comparison the left plot shows the con^esponding prediction at NLO accuracy using the PDF set CTEQ6.5 
[24]. 



4 Conclusions 

In this article we have summarized the present knowledge on theory predictions for the top-quark 
pair production cross section at Tevatron and LHC. We have taken some care to quantify the sen- 
sitivity of the total cross section to soft gluon emission and large Sudakov type logarithms. As 
is well known, top-quark pair production at Tevatron is largely dominated by parton kinematics 
close to threshold, thus approximations based on soft gluon resummation should provide an ex- 
cellent description. At LHC we find that soft gluon emission near threshold is less dominant, but 
contributes still a numerically sizable fraction to the total cross section. Thus, soft gluon effects in 
/^/"-production are still rather prominent at LHC as well. 

We have updated the NLL resummed cross section as defined in Ref. [8, 14] using recent PDFs. 
Furthermore, we have extended the resummed predictions to NNLL accuracy and we have derived 
approximate NNLO cross sections which are exact to all powers in ln|3 at two loops. Together 
with the exact NNLO scale dependence (and including the two-loop Coulomb corrections) our 
result for OjsjjvfLo (approx) represents the best present estimate for hadro-production of top-quark 
pairs, both at Tevatron and LHC. As mentioned earlier we believe that hard corrections at the 
NNLO level are small. This is supported by the explicit findings of Ref. [15]. We have found 
for the NNLL resummed cross section and the finite order expansion good apparent convergence 
properties. iVloreover, the stability of the total cross section with respect to scale variations is much 
improved by our NNLO (approx) result. 

In closing let us briefly comment on ideas to use top-quark pair production as an additional 
calibration process for the parton luminosity at LHC [11]. This could become feasible because 
the PDF dependence of ff-production at LHC is anti-correlated with W /Z-hoson production (the 
standard candle process at LHC, see e.g. [42,43]) and correlated with Higgs boson production, 
especially for larger Higgs masses. It has been noted however, that the NLO theory predictions 
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to the top-quark cross section are not accurate enough. We are confident that the NNLO (approx) 
resuks of this present paper provide a step in the right direction by further constraining the theory 
uncertainties for this important process. 

The complete NNLO QCD predictions for heavy-quark hadro-production do not only require 
the hard scattering cross section but also the evolution of the parton densities to be performed at 
the same order employing the NNLO splitting functions [29, 30]. In addition we remark, that the 
present accuracy on the gluon PDF in the medium x range of interest for top-quark production at 
LHC is well constrained by from DIS data for structure functions in e"*^/?- scattering from HERA 
and evolution, leading to the rather small uncertainty of 3% (see Fig. [I]) at small energies. Thus, 
for top-quark pair production to become a standard candle process at LHC similar to W /Z gauge 
boson production and to become competitive with DIS data also an experimental accuracy much 
better than the currently quoted [1] value of 10% will be needed. 
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A Useful formulae 

The Bom cross section in the color-basis defined by color-singlet and color-octet final states reads 
in momentum space with p = Am^ j s and P = ^1 — p 
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where Ca and Cf are the usual color factors, with C^= N = 3 and Cf = jiq{N^ ~ 1) = 4/3 in 
QCD. The Mellin moments as defined in Eq. ([7]) can be easily computed, see e.g. Ref. [8]. 

Next we present the perturbative expansions for the anomalous dimensions Aq, Ag, Dq and Dg 
entering Eq. (fT3l) . We have for the quark case [29, 37] 
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where nf denotes the number of effectively massless quark flavors and ^2, ^3, • • ■ are the values of 
the Riemman zeta-function. Likewise, the Dq read 
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All gluonic quantities are given by the simple relation 

Ai) _ Ca (i) 
^8 ~ r 1 ' 
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Here we summarize the functions gjj, gj, j, g], j appearing in the resummed cross section 
Eq.lfTOl) to NNLL accuracy [26-28]. Keeping the full dependence on fj^ and fir, we have 
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with Ye = 0.5772167. The gluonic expressions gip,gl„ i and g^ are obtained with the obvious 

replacements Aq -^ Ag and Dq — »• D^ . The dependence on |3o is recovered by A*-'-* -^ A^'7 |3q, 
d(') -^ D*^')/(3q, P/ ^ (3,7(3q^^ and multiplication of g]: j by (3o. We also give explicit results for 
the matching functions gj* j in Eq.®, 
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where C^J^ = 367ra° + (n^ - 4) (2/31n2 - 5/9) and C^J^ = 768/77la° with the numerical constants 
aj = 0.180899 and a° = 0.108068 being reported in Table 1 of Ref. [4]. The presently unknown 

(2) (2) 

two-loop constants are denoted Q=^ and Cgg . 

At first and second order in a^ the Coulomb corrections have to be added to the cross section. 
In the limit p ^ 1 (that is P ^ 0) they read [33, 38] for the color-singlet final state 
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From these expressions the octet results for a)- g and 6^ ■ g are obtained by the replacement of 
the color factor 2Cp -^ (2C^ ~C^) consistent with results from potential non-relativistic QCD 
(e.g. [38,39]). 

Next we give the Mellin transforms of powers of logarithms in [3 according to Eq. (O. Recall 
that p = Amt'^ /s and |3 = \/\ —Amp-js = y/1 —p. In the limit p ^ 1 as needed in the fixed order 
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where A'' = A''exp(Ye). 
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and Cgg has again been given below Eq. (IAI13I) while Qg is the yet uncalculated two-loop con 



stant. 
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Table 10: The total cross section for /j = mt at LHC for nit = 171 GeV and the full set of predictions 
from the MRST-2006 NNLO PDF set [25]. All rates aie in pb. We denote by a^LQ [4, 5] the NLO QCD 
prediction by ares the result of NLL threshold resummation [8] and by <7nnlo (approx) ^^^ NNLO QCD 
prediction based on soft gluon approximation and exact two-loop scale dependence [9]. The upper and 
lower indices denote the shifts towards /j = 2m, and ^ = m,/2. 
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Table 1 1 : The total cross section for /j = nit at LHC for m, = 171 GeV and the full set of predictions from 
the CTEQ6.5 PDF set [24]. All rates are in pb. We denote by a^LQ [4,5] the NLO QCD prediction by ares 
the result of NLL threshold resummation [8] and by CJjyjj^LQ (approx) ^^ NNLO QCD prediction based on 
soft gluon approximation and exact two-loop scale dependence [9]. The upper and lower indices denote the 
shifts towards /u = 2m, and /u = m,/2. 
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Table 12: The total cross section for ^ = m, at Tevatron for nit = 171 GeV and the full set of predictions 
from the CTEQ6.5 PDF set [24]. All rates are in pb. We denote by a^LQ [4,5] the NLO QCD prediction 
by CJres the result of NLL threshold resummation [8] and by Oj^nlq (approx) ^^ NNLO QCD prediction 
based on soft gluon approximation and exact two-loop scale dependence [9]. The upper and lower indices 
denote the shifts towards /j = 2m, and ^ = m,/2. 
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30 


7 64-"'^'' 

'■0^+0.47 


^■83;°:- 


Q 9Q-0.34 
^■2«+0.08 



Table 13: The total cross section for ^ = mt at Tevatron for nit = 171 GeV and the full set of predictions 
from the MRST-2006 NNLO PDF set [25]. All rates are in pb. We denote by a^LQ [4, 5] the NLO QCD 
prediction by ares the result of NLL threshold resummation [8] and by CT|yfjyfLQ (approx) ^^^ NNLO QCD 
prediction based on soft gluon approximation and exact two-loop scale dependence [9]. The upper and 
lower indices denote the shifts towards /u = 2m, and /a = m,/2. 
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